Emergent patterns in a spin-orbit coupled spin-2 Bose-Einstein condensate 
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The ground-state phases of a spin-orbit (SO) coupled atomic spin-2 Bose-Einstein condensate 
(BEC) are studied. Interesting density patterns spontaneously formed are widespread due to the 
competition between SO coupling and spin-dependent interactions like in a SO coupled spin-1 con- 
densate. Unlike the case of spin-1 condensates, which are characterized by either ferromagnetic 
or polar phase in the absence of SO, spin-2 condensates can take a cyclic phase, where we find 
the patterns formed due to SO are square or triangular in their spin component densities for axial 
symmetric SO interaction. Both patterns are found to continuously evolve into striped forms with 
increased asymmetry of the SO coupling. 
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I. INTRODUCTION 



The technique of gauge field is a useful tool in theo- 
retical physics. Depending on commutation relations of 
the associated operators, gauge fields are classified into 
Abelian and non-Abelian ones. An important example 
for Abelian gauge field concerns the well studied vector 
and scalar potentials of electromagnetic fields. Interest- 
ing phenomena, such as integer and fractional quantum 
hall effects are observed in two-dimensional high mobility 
electrons inside a perpendicular magnetic field. For a col- 
lection of cold atoms, Abelian gauge field can be affected 
through rotation [J B or adiabatic motion inside far off- 
resonant laser fields [3|-[6[ . While Abelian gauge fields are 
often studied, recent interests in Bose condensed atoms 
are increasingly targeted at situations when gauge fields 
become non-Abelian [H, 0-@- One of the simplest ex- 
amples concerns spin-orbit (SO) coupling, which recent 
understandings increa singly view as exists intrinsically 
in solid-state systems [ltj, El i and is responsible for the 
quantum spin hall effect. In atomic quantum gases, SO 
coupling has been proposed with a variety of approaches 
involving atomic interactions with electromagnetic fields 
0, IT2h18| . including schemes for two component or 
pseudo spin- 1/2 atomic condensates [13, [lj| and spin- 
1 condensates @, HH, and realized experimentally for 
pseudo spin- 1/2 system recently [lq ]. 

Spin-orbit interaction couples the internal (spin) and 
orbital (momentum or angular momentum) degrees of 
freedom. In cold atoms, elastic binary collisions are de- 
scribed by contact interactions proportional to their re- 
spective s-wave scattering lengths. For atoms with inter- 
nal degrees of freedom, their collision interactions gen- 
erally contain spin dependence, as in spinor atom con- 
densates [l9| ] . where the spin-dependent interactions de- 
termine the various ground state phases: ferromagnetic 
or polar for the case of spin-1; ferromagnetic, polar, or 
cyclic, for spin-2; etc. The inclusion of SO coupling in- 
duces competition with spin-dependent interactions in 
addition to modify the single particle spectra for a spinor 
BEC. As a result, spin-dependent interaction will in turn 



influence atomic spatial motion, leading to a variety of 
density patterns even in the ground state fl3l - [l6| . 

Spontaneously formed patterns in spinor component 
densities, are essentially spin textures, previously stud- 
ied in 3 He superfluids and more recently in spinor con- 
densates for both ferromagnetic and anti-ferromagnetic 
(polar) phases, especially with long-range dipolar inter- 
actions |20l-|26j. Periodically ordered patterns can also 
arise in superfluids with roton-like spectra [H, [24| • 

This study addresses ground state pattern formations 
due to the competition between SO coupling and spin- 
dependent interactions. The simplest case concerns a 
spin-1/2 BEC with SO coupling, Qj|{l5|, already ob- 
served in recent experiments [16[ . Depending on the sign 
of the effective spin interaction strength, which is propor- 
tional to the difference between intra- and inter-species 
atomic scattering strengths, the ground state patterns 
are found to be either planar or standing waves EMI. 
The same conclusion was reached in recent theoretical 



studies for a spin-1 BEC (13|, where the two ground state 
phases are ferromagnetic and polar in the absence of SO 
coupling. Spin-2 condensates, on the other hand, can po- 
tentially be very different due to the presence of a cyclic 
phase when SO coupling is absent. Understanding the 
associated patterns in a spin-2 BEC thus constitutes an 
importance objective for new physics in atomic quantum 
gases with gauge fields. 



II. OUR MODEL 

We start by describing our model of a spin-2 BEC 
[27l |28| in a quasi-two dimensional optical trap V D = 
mw 2 (x 2 + y 2 + \ 2 z 2 )/2 with A > 1, including SO cou- 
pling. The effective Hamiltonian takes the form 



H = J dpi\}\ + V x p x F x + VyPyFy + t/j J 

+ 1 '(-iy+^UU^j}: (i) 
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where {F^)ij (fi — x,y,z) are the components of 

spin-2 matrices F^, V' Q = mu 2 (i 2 + y 2 )/2, a', /?' and 
7' are reduced two-dimensional effective density-density, 
spin-exchange, and spin singlet-pairing interaction pa- 
rameters, respectively. v x and v y parameterizes the 
strength of SO coupling. This form of two- and three- 
component of condensates with SO coupling can be real- 
ized with a tripod and a tetrapod plane- wave laser beams 
setup 0], respectively. Also, recently Y. - J. Lin et al. [HI 
have realized a mixed form of Rashba and Dresselhaus 
type SO coupling interaction for pseudo spin-1/2 sys- 
tem. They further suggested that we can realize the SO 
coupled spin-1 or spin-2 condensates by using a smaller 
quadratic Zeeman shift. 

The single particle Hamiltonian 

P 2 

HQ = + V X p X F X + VyPyFy, , (2) 

ignoring the external trapping potential, can be easily di- 
agonalizcd using plane wave basis, leading to eigenvalues 
E\(p) and eigenvectors <j)\ (k) = e p £\(k) given by 

^A(k) = ^+A V /^fc2 +rA ,2 fc 2 5 (3) 

C 2 T (k) = (xk 2 ,2 X £,V6,2 Xk , x 2 )/4; 

d T (k) = (-XkVxk,o,x k ,Xk)A 

C T (k) - (^ Xk 2 ,0,-\/2,0,^Xk)/2\/2, 

C T !(k) = (-Xk 2 ,Xk,o,-xk,x k )A 

C- 2 (k) = (XkV2x k ,%/6,-2xk,Xk)/4. (4) 

In the above A = ±2, ±1,0 labels the respective energy 
band and Xk = {v x k x + iv y k y ) / 'v x k x + v%k%. Cor- 
responding to the above eigenvectors Eq. Q, we find 
the two order parameters (|(F)| = |C^FCa| and |(6)| = 
\{-VCx.jCx,-j)\) are equal to (2,0), (1,0), and (0, 1) cor- 
responding to A = ±2, ±1, 0, respectively. 

The ground state of the single particle Hamiltonian 
([2]): </>_2(k s ), is two-fold degenerate for \v x \ ^ \v y \; and 
is infinitely degenerate if \v x \ = \v y \. For the former 
case, hkg = (±2m\v x \,0) if \v x \ > \v y \, while Kk g = 
2m\v x \{cos9 g , smOg) for any 9 g <G [0, 2n), in the latter 
case with \v x \ = \v y \. Atomic spins are fully polarized 
with |(F) | = 2 for both cases. 



it is weakened apprcciabely for numerical studies with 
imaginary-time propagation of coupled Gross-Pitaevskii 
equations (GPEs). Including axisymmetric SO interac- 
tion, the ground states can be approximately understood 
by applying perturbation theory to the degenerate single 
particle states. Each of the three types of spin-dependent 
interactions is found to be associated with a distinct num- 
ber single-particle states. For the case of "cyclic" inter- 
action, this approximation breaks with increasing asym- 
metry of SO coupling. 

(A) "Ferromagnetic" interaction The ferromagnetic 
spin-exchange interaction tends to polarize atoms' spin, 
consistent with the ground state of the single-particle 
Hamiltonian. As a result, all atoms condense into any 
single particle ground state 0_2(k s ). The spin rotation 
and time-reversal symmetry are broken in this case. The 
optical trap and the spin-independent interactions will 
change atomic density distributions slightly. 




FIG. 1: (Color online). Ground states of a spin-2 BEC in 
the "polar" interaction case with SO coupling for v x = 2v y 
and 7' = — P' = —0.2a' < 0. (a-e) Density distributions for 
the five spin components: Mf = +2, +1, 0, — 1, — 2. Yellow 
and black colors represent high and low density respectively, 
(f-g) Spatial distributions of two order parameters (uniform 
in this case) |(F)| and |(G}|, which approximately measure 
spin-exchange and singlet-pairing interactions, respectively. 



III. GROUND-STATE PATTERNS 

When SO coupling is absent, the ground states for a 
spin-2 BEC including collision interactions, takes three 
phases: ferromagnetic, polar, and cyclic [23, Their 
corresponding spin-dependent interactions will in this 
study be called "ferromagnetic", "polar", and "cyclic" 
interactions. In the presence of SO coupling, we dis- 
cuss their possible ground states for the above three 
types of spin-dependent interactions. The external trap- 
ping potential is neglected in analytical treatments, while 



(B) "Polar" interaction The order parameter 
|(F) I = 2 in the single particle ground state conflicts 
with the polar phase of |(F)| = 0. A superposition of 
degenerate single particle states is adopted 

V.> = ]Tc g ^ 2 (k ff ), (5) 

to minimize the spin-dependent interaction. The states 
</>_2 (k g ) arc orthogonal, thus the superposition state 
docs not affect kinetic energy. The suitable coeffi- 
cients C g should give \{ip\F\i())\ = and K^l©!^)! = 1, 
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FIG. 2: (Color online). Ground states of a spin-2 BEC in the 
"cyclic" interaction case with SO coupling for v x = v y . (a) 
The number of plane waves in the ground states with changing 
if = arg(7' + i/3'). The column of 4 blue (3 red) dots denote 
four (three) plane waves. The shadow window with one or two 
black dots denote one or two plan waves respectively consti- 
tuting of the ground states. Schematic illustrations of four 
b) [three c)] plane wave superpositions responsible for square 
(triangular) density modulations. The angles between the di- 
rections of two nearby plane waves are close but not exactly 
equal to 27r/4 or 2tt/3 due to frustration. 



while preserving the spin-independent interaction energy. 
It is easy to conclude a superposition of two counter- 
propagating single particle states meets this criterion, 



thus the ground state is given by 0_2(k s ) + e 



-k s ), 



where the time-reversal symmetry remains valid. Back 
into the z-axis quantized representation, the interference 
between the two counter-propagating plane waves induce 
density oscillations in spinor components. In the quanti- 
zation direction of k ff , however, only Mp = ±2 compo- 
nents are found to be populated, each with smooth den- 
sity, respectively propagating into opposite directions. 

In Figure [1] we illustrate density distributions of the 
spin components and the spatial dependence of two order 
parameters |(F)| and |(9)| for v x = 2v y and 7' = — (3' = 
-0.2a' < 0. When \v x \ > \v y \ (\v x \ < \v y \), we find 
the plane waves propagate along the positive/minus x- 
axis (y-axis) the ground state, resulting in density mod- 
ulations along x-axis (y-axis). Figure [TJf-g) shows the 
two order parameters are uniformly distributed, with 
|(F) I = and |(0)| = 1, consistent with previous dis- 
cussions. A cautionary note concerns the axis-symmetric 
case of \v x \ = \v y \. Due to the spin-independent interac- 
tion, the ground state consists of only one pair of counter- 
propagating plane waves, although superpositions of two 
or more counter-propagating pairs also minimize spin- 
dependent energy. 

(C) "Cyclic' interaction First, we consider the sym- 
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FIG. 3: (Color online). Ground states of a spin-2 BEC in 
the "cyclic" interaction case with SO coupling for v x = v y , 
7' = 0.2a', and j3' — 0.02a'. (a-e) Density distributions of the 
spin components: Mf = +2, +1,0, —1, —2. Yellow and black 
colors represent high and low densities respectively, (f-g) Spa- 
tial distributions of the order parameters |(F)| and |(6)|, re- 
flecting spin-exchange and singlet-pairing interactions respec- 
tively. 
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FIG. 4: (Color online). The same as in Fig. [3]but for : 
7' = 0.02a' and /3' = 0.2a'. 



metric case of \v x \ = \v y \. As pointed out in the above 
"polar" interaction case, superpositions of single par- 
ticle ground states are tried to minimize the interac- 
tion energy. Numerically, the minimization procedure 
is achieved with simulated annealing. In Fig. EJa), 
the ground states are found to contain several plane 
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FIG. 5: (Color online). Ground states of a spin- 
2 BEC in the "cyclic" interaction case with SO 
coupling for 7' = = 0.2a'. (al-a5) Spinor com- 
ponent densities at v y — v x , from up to down 
corresponding to Mf — +2, +1,0, 1,2. (bl-b5), 
(cl-c5), (dl-d5), and (el-e5) are as in (al-a5) but 
for v y = 3v x /4, v y = 17^/24, v y = 2v x /3, and 
Vy = v x /2, respectively. 



waves with significant weights if they are assumed ap- 
proximately to consist of a series of degenerate sin- 
gle particle ground states with momentum hk. g = 
2m\v x \(cos9g,sm.9g) and if external optical trapping po- 
tentials are neglected. In this case, the spin-dependent 
interactions couple degenerate single particle ground 
states, and result in discrete rotation symmetry of the 
ground states. Starting from the boundary of /?' = 0, 
with increasing value of tp = arg(7' + i/3'), the ground 
states at first contains only one plan wave, as tp ex- 
ceeds a small critical value determined by the three in- 
teraction parameters, the ground states display periodic 
square lattices in the spin component densities. Further 
increasing of tp arrives at the critical value tp c , where the 
ground states change into triangular from superpositions 
of three plane waves until near the boundary with 7' = 0, 
where the ground state superpositions reduce to contain 
two counter-propagating plan waves. It is important to 
emphasize that interference between serval plane waves 
may increase density-density interaction energy, which 
explains why the phase boundaries between the ground 
states with four and one, also with three and two plan 
waves, deviate from tp = or ip = n/2. The numerical 
results show for the two different types of ground states 
with square or triangular patterns, the angles between 
nearby plane waves are close to but not exactly equal to 
27r/4 or 27r/3 respectively, and the nonzero coefficients C g 
are not all exactly equal, although they are of the same 
orders of magnitudes. 

To validate our proposition that the ground states of 



a spin-2 BEC with SO interaction are approximate su- 
perpositions of degenerate single particle ground states 
in the absence of the trap potential, we numerically solve 
the GPEs directly using imaginary-time propagation in- 
cluding the optical trap. Figures [3] and 0] show spin com- 
ponent densities and the spatial dependence of the two 
order parameters | (F) | and | (0) | for two different cases. 
They confirm our understanding that not only spin com- 
ponent density distributions but the two order parame- 
ters show square or triangular patterns. More carefully, 
the shapes of density distributions within each unit cell 
are not identical for the five spin components at a fixed 
tp. Their differences vary as tp is changed. Because the 
values of C g are not exactly equal, they generally vary 
with tp constrained by the spin-independent interaction 
term, although the number of plane waves contained in 
the ground states remains at four or three. 

Next, we consider the asymmetric case of \v x \ 7^ \v y , 
where the single particle ground states are two-fold de- 
generate. When insufficient number of degenerate single 
particle ground states exist to construct square or trian- 
gular shaped density distributions, our proposition fails. 
This is also validated by numerical simulations. In Fig. 
we illustrate results from numerical simulations. From 
(al-a5), (bl-b5), (cl-c5), (dl-d5), to (el-e5), the asym- 
metry of the SO coupling is increased, realized by fix- 
ing v x while decreasing v y from being equal v x , 3^/4, 
171^/24, 2-1^/3 to v x /2, respectively. The ground state 
show triangular density distributions in each component 
at v y = v x . The triangular lattice is deformed into al- 
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most square lattice distributions when v y is decreased 
into 3v x /4, and gradually evolve in the end into striped 
patterns reflecting the underline two-fold degeneracy as 
we further decrease v v . 



IV. CONCLUSION 

In conclusion, we study spin-2 condensates with SO 
coupling. Due to the competition among SO coupling, 
spin-dependent, and spin-independent interactions, the 
ground states are found to contain one or two counter- 
propagating plane waves respectively in "ferromagnetic" 
and "polar" interaction cases, where both the single par- 
ticle Hamiltonian and atomic interaction energies can 
reach their corresponding minimum values. For "cyclic" 
interaction case with axisymmetric SO coupling, the 
ground states can be categorized into two different types 
respectively containing four or three plane waves, lead- 
ing to square or triangular patterns correspondingly. For 



asymmetric SO coupling, the ground states are deter- 
mined by the ratio of two SO coefficients: |%|/|wa:- 

As 

long as ratio is decreased, the spin component density dis- 
tributions can evolve from triangular to square shaped 
lattices before asymptotically reaching the striped pat- 
tern. This structure phase transition provides a clear 
signature for SO coupling. In the "cyclic" interaction 
case, the corresponding minimum value of single parti- 
cle Hamiltonian and the interaction energies cannot be 
reached simultaneously. 
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